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Abstract — The multicast capacity of a directed networls is 
closely related to the s-t maximum flow, which is equal to 
the s-t minimum cut capacity due to the max-flow min-cut 
theorem. If the topology of a network (or link capacities) is 
dynamically changing or have stochastic nature, it is not so trivial 
to predict statistical properties on the maximum flow. In this 
paper, we present a coding theoretic approach for evaluating 
the accumulate distribution of the minimum cut capacity of 
weighted random graphs. The main feature of our approach 
is to utilize the correspondence between the cut space of a graph 
and a binary LDGM (low-density generator-matrix) code with 
column weight 2. The graph ensemble treated in the paper is a 
weighted version of Erdos-Renyi random graph ensemble. The 
main contribution of our work is a combinatorial lower bound 
for the accumulate distribution of the minimum cut capacity. 
From some computer experiments, it is observed that the lower 
bound derived here reflects the actual statistical behavior of the 
minimum cut capacity. 

I. Introduction 

Rapid growth of information flow over a network such as 
a backbone network for mobile terminals requires efficient 
utilization of full potential of the network. In a multicast 
communication scenario, it is well known that an appropriate 
network coding achieves its multicast capacity. Emergence of 
the network coding have broaden network design strategies for 
efficient use of wired and wireless networks [16]. 

The multicast capacity of a directed graph is closely related 
to the s-t maximum flow, which is equal to the s-t minimum 
cut capacity due to the max-flow min-cut theorem [15]. The 
topology of a network and the assignment of the link capacities 
determine the s-t minimum cut capacity of a network. 

If a network is fixed, the corresponding s-t maximum 
flow of the network can be efficiently evaluated by using 
Ford-Fulkerson algorithm [15]. However, if the topology of 
a network (or link capacities) is dynamically changing or 
have stochastic nature, it is not so trivial to predict statistical 
properties on the maximum flow. For example, in a case of 
wireless network, the link capacities may fluctuate because 
of the effect of time-varying fading. Another example is 
an ad-hoc network whose Unk connections are stochastically 
determined. 

In order to obtain an insight for the statistical property of 
the min-cut capacity for such random networks, it is natural to 
investigate the statistical properties of a random graph ensem- 
ble. Such a result may unveil typical behaviors of minimiun 



cut capacity (or maximum flow) for given parameters, such as 
the number of vertices and edges. 

Several theoretical works on the maximum flow of random 
graphs (i.e., graph ensembles) have been made. In a context of 
randomized algorithms, Karger showed a sharp concentration 
result for maximum flow in the asymptotic regime [11]. 
Ramamoorthy et al. presented another concentration result; 
the network coding capacities of weighted random graphs and 
weighted random geometric graphs concentrate around the 
expected number of nearest neighbors of the source and the 
sinks [12]. These concentration results indicate an asymptotic 
property of the maximum flow of random networks. Wang 
et al. shows statistical property of the maximum flow in an 
asymptotic setting as weU. They discussed the random graph 
with BemouUi distributed weights [10]. 

In this paper, we present a coding theoretic approach for 
evaluating the accumulate distribution of the minimum-cut 
capacity of weighted random graphs. This approach is totally 
different from those used in the conventional works. The 
basis of the analysis is the correspondence between the cut 
space of an undirected graph [5] and a binary LDGM (low- 
density generator-matrix) code with colunm weight 2. Yano 
and Wadayama presented that an ensemble analysis for a class 
of binary LDGM codes with column weight 2 for the network 
reUability problem [13]. This paper extends the idea in [13] to 
weighted graph ensembles. We focus on a weighted version 
of Erdos-Renyi random graph ensemble [4] in this paper. 

II. Preliminaries 

In this section, we first introduce several basic definitions 
and notation used throughout the paper. Then, the cut-set 
weight distribution wiU be discussed. 

A. Notation and definitions 

A graph G = {V. E) is a pair of a vertex set V = 
{wi, . . . , Vk} and an edge set E = {ei, . . . , e„} where Cj = 
{u,v),u,v e y is an edge. If ej = {u,v) is not an ordered 
pair, i.e., {u,v) = {v,u), the graph G is called an undirected 
graph. Otherwise, i.e., [u, v) is an ordered pair, G is a directed 
graph. The two vertices connecting an edge e G E are referred 
to as the end points of e. If an edge e = {u, u) has the identical 
end points, e is called a self-loop. 



If real valued function w : E R>o is defined for 
an undirected graph G = {V,E), the triple {V,E,w) is 
considered as a weighted graph. The set M>o represents the 
set of non-negative real numbers. In our context, the weight 
function w represents the Unk capacity for each edge. 

Assume that a weighted undirected graph G = {V,E,w) 
is given. A non-overlapping bi-partition V = X U {V\X) is 
called a cut where X is a non-empty proper subset of V{X ^ 
y). The set of edges bridging X and V\X is referred to as the 
cut-set corresponding to the cut of {X, V\X^. The cut weight 
of X is defined as E„ex,«ey\x ^("' 

B. Random graph ensemble 

In the following, we will define an ensemble of weighted 
undirected graphs. The graph ensemble is based on Erdos- 
Renyi random graph ensemble. Let k(k > 1) be the number 
of labeled vertices and n(l < n < k(k — l)/2) be the number 
of labeled undirected edges. The vertices are labeled from 1 
to k and the edges are labeled from 1 to n. 

For any adjacent vertices, a single edge is only allowed. It 
is assumed that each edge has own integer weight; namely, 
a weight Wi E [i-.q](i G [1/"]) is assigned to the edge with 
label i, which is denoted by the ith edge. The notation [a, b] 
denotes the set of consecutive integers from a to b. The set 
i?^ ^ denotes the set of all undirected weighted graphs with 
fc-vertices and n-edges satisfying the above assumption. 

For any G G Rl^, the sets of vertices and edges are denoted 
by V{G) and E{G), respectively. In a sinoilar way, Wi{G) is 
defined as the weight of ith edge of G. 

It is evident that the cardinality of i?^ ^ is given by 



\Rln\ 



(2) 



(1) 



We here assign the probabiUty 

P{G) = — ^mK(G))/xK(G)) • • • /xK(G)) (2) 

for G e where ^ is a discrete probability measure defined 
over [l,q]; namely, it satisfies 

^ /x(u;) = 1 and Vw e [1, q\,IJ.{w) > 0. (3) 

The pair (i?^ n^^) defines an ensemble of random graphs 
and it is denoted by £. 

C. Incidence matrix 

For G G Rkrv '^he incidence matrix of G, denoted by 
M(G) e {0, Ijfc'x", is defined as follows: 



M{G\i 



1, if i-th vertex connects to j-th edge 
0, otherwise, 



(4) 

where M{G)ij is the (z, j)-element of M{G). If G is con- 
nected, then the rank (over F2) of M{G) is — 1. The row 
space (over F2) of M{G) coincides with the set of all possible 
incidence vectors of cut- sets of G. 



D. Cut weight distribution 

For a given undirected graph, we can enumerate the number 
of cut- sets with cut weight w. The cut weight distribution of 
Gby 

Bw{G) = ![£' is a cut-set of G,cut weight is w] 

E<ZE(G) 

(5) 

for positive integer w. The function ![•] is the indicator 
function that takes value 1 if the condition is true; otherwise it 
takes value 0. This cut weight distribution can be regarded as 
an analog of the weight distribution of the binary linear code 
defined by the incidence matrix of a given undirected graph. 

For ensemble analysis, it is convenient to introduce another 
form of the weight distribution. The detailed cut weight 
distribution Au,v.w '■ Rtn^ ^>o is defined by 



E E 



iM{G) 



c,Y^CiWi{G) = w 

1=1 



,(6) 



for M G [1 , fc — 1] , w G [0, n] . The set of constant weight binary 
vectors Z^" ''^ is defined as 



Z^'^''') = {xG {0,1}" ■.WH{x)^b}. 



(7) 



The function wh{-) represents the Hamming weight. Assume 
that the cardinality of the cut is v and that the size of X c .B 
is u. Under this condition, the function Au,v.w{G) represents 
the number of cuts with the cut weight w. It should be noted 
that the one-to-one correspondence between the cut space and 
the set of incident vectors of the cuts are impUcitly used in 
the definition of Au,v,w{G). 

The following lemma indicates the relationship between 
and B^{G). 

Lemma 1: For G G i?^ the cut weight distribution 
Bw {G) can be upper bounded by 

k—l n 

(G), (8) 

for w e Z>o. The notation Z>o represents the set of non- 
negative integers. 
Proof: Let 

S{G) = {c = mM{G) G I m G F^\{0^ l'^}}. 
The cut weight distribution B^iG) can be rewritten as follows. 
Bw{G) = is a cut-set of G, cut capacity is ti;] 

E(1E(G) 



E ^ 

ceS(G) 



YciWi{G) = w 

Li=l 



(9) 



The second equality is due to the fact that the row space of 
M (G) equals the set of all possible cut-set vectors of G. It is 
evident that 

|{m G F^ I TO 9^ 0*^, m ^ 1*^, c = toM(G)}| > 2 



holds for any c £ S{G). This impHes that 



M (G). Since M{G) is an incidence matrix, the column weight 
of fi is WHifi) — 2 for i £ [l,n]. From the assumptions 
^ h{c)l [c = mM(G)] described above, we have 



ceS(G) 



me{o,i}''\{o'=,i'=} ce{o,i}" 



(10) 

holds for any real-valued function h : {0, 1}" — > M. Substitut- 
ing ( [T0| into (|9]l, we obtain 

BUG) < I J2 E I[c-"^M(G)] 



me{o,i}''\{o'=,i''} ce{04}" 



X I 



,i=l 
_j ,c— 1 n 

oEE^«.-(G) 



(11) 



11=1 ^=0 



III. Ensemble average of cut weight distribution 

In this section, we discuss the average of Au.v,w{G) over the 
ensemble £. This analysis is very similar to the derivation of 
the average weight distribution of LDGM codes with column 
weight 2. 

A. Preparation 

In the following, the expectation operator E is defined as 

E[/(G)]- J2 (12) 

where / is any real-valued function defined on Rj, ^. The next 
lemma plays a key role to derive a closed form the average 
cut set weight distribution. 

Lemma 2: Assume that m* £ {0, 1}*^ and c* £ {0, 1}" 
satisfies WHim*) — u and wh{c*) — v where u £ [1, fc — 1] 
and V £ [0,n]. The following equality 



' M{G) = c* ,J2'^>^iG) ^ w 



u{k-u)\ fCl) -u{k-u) 



0{% 



holds. The function f{x) is defined by 



(13) 
(14) 



The term [x'^]f{x)'" represents the coefficient of x'" in f{x)'". 

Proof: Due to the symmetry of the ensemble, we can 
assume that the first w-elements of m* are one and the rests 
are zero without loss of generality. In a similar manner, c* is 
assumed to be the binary vector such that first w-elements are 
one and the rests are zero. 

In the following, we will count the number of labeled 
graphs satisfying m*M{G) = c* by counting the number of 
binary incidence matrices satisfying the above condition. Let 
M{G) = (/i /2 ■ ■ ■ fn) where is the ith column vector of 



"^*/. = ( n '^T'i 1 (15) 
\ 0, ^ e [u + l,n]. ^ ' 

We then count the number of allowable combinations of 



(./i, /2, • • ■ , fn) satisfying ( 15 i. Let 



A^{f £{QAf\m*f ^l,WH{f)^2). 



(16) 



The cardinality of A is given by ~ u{k—u) because a non- 
zero component of / needs to have an index within [l,u] and 
another non-zero component has an index in the range [u + 
l,k]. This observation leads to the number of possibilities for 
ifi, f2, ■ ■ ■ ; fv) which is given by The remaining 

n — V columns, {fv+i, ■ ■ ■ , /n), should be taken from the set 
{/ £ {0, l}*^ I wnlf) = 2}\^. Thus, the number of possibil- 
ities for such choice is (n — v) . In summary, the 
number of allowable combinations of (/i, /2, . . . , /„) denoted 
by S is given by 



S — vl(n — w)! 



u{k — u)\ / (2) — u{k — u) 



(17) 



We are now ready to derive the claim of this lemma. To 
simplify the notation, the cut weight is denoted by tp(c) — 
J27=i Ci'^i{G)- The left hand side of (13 1 can be rewritten as 
follows: 

E[I[m*A/(G) =c*,V(c*) = 

= P{G)I[m*M{G) = c*Xc*) = w] 

1 



n 



.n, E n 

Vn) GeRl „ \^ell,n] 



X I[m*M{G) = c*,'iJj{c*)^w] 



E 




«i+2m2H — quq=w 

)\(u{k-u)\({'^^-u{k~u) 



n — v 



[x^]f{x)\ 
(18) 



The last equality is due to ( [T7] l. ■ 

The following lemma provides the ensemble average 
^[Au,v,w{G)\, which is a natural consequence of Lemma |2] 
Lemma 3: The expectation of Au^v,w{G) is given by 



EK,„,^(G)] 

^ ((2)) W V V )\ n~v 
where u £ [l,k — l\,v £ [0, n], w £ Z>o- 



[x^]f{x)\ (19) 



Proof: The expectation of A„ „^(G) can be simplified 
as follows: 



= E 



^ ^ l[mM{G) = c, 0(c) = w] 
^ ^ £[l[mM{G) = c,4,{c)=w]] 

'"^E[I[m*M(G) = c*,VXc*) (20) 



The last equality is due to the symmetry of the ensemble. The 
binary vectors m* € {0, 1}'' and c* G {0, 1}" are arbitrary 
vectors satisfying wni'm*) = u and wjy(c*) = v. Substituting 



( [T3) l in the previous Lemma into (20 1, we obtain the claim of 
this lemma. ■ 

B. Upper bound on average cut weight distribution 

In order to investigate statistical properties of the minimum 
cut weight, it is natural to study the tail of the average cut 
weight distribution. The following theorem provides an upper 
bound on average cut weight distribution that is the basis of 
our analysis. 

Theorem 1: The expectation of Bw{G) over £ can be upper 
bounded by 

2(y) ^1 v^; „t^o V V 
'(2) ~ "^(^ ^ 



for w e Z>o. 

Proof: Due to Lemma [T[ we immediately have 

k-l 

E[B^(G)] < E ' ' 



(21) 



2 ^ ^ Au^v.w{G) 

u—l v—Q 
-. n,— 1 n 

„^^EK,_(G)]. 



(22) 



(23) 



u—l v—0 



Substituting the left hand side of the equality in Lemma [3]into 
the last equation gives the claim of the theorem. ■ 

C. Accumulate cut distribution 

Let us define the accumulate cut weight of G, Gs{G), by 

<5-l 

Cs{G) - ^ BUG), (24) 

to=0 

where J is a non-negative integer If Cs (G) is zero, the graph 
G does not contain a cut with weight smaller than 6. This 
implies that A(G) > 5 in such a case, thus we have 



Pr[X{G) > S] 



Pr[C5iG)^0] 
l-Pr[Cs{G) > 1]. 



The second equality is due to the non-negativity of Cs{G). The 
probability Pr[X{G) > 6] can be considered as the accumulate 
probability distribution for the minimum cut capacity: 



Pr[\{G) >6]= J2 ^(G)I[A(G) > S] 



(26) 



GERl 



The following theorem is the main contribution of this work. 

Theorem 2: Assume that an ensemble £ is given. The 
probability Pr[X{G) > 6] can be lower bounded by 



-. (5— 1 fc — 1 / 1 \ n 



w—O u—l 



v=0 



u{k — u) 

V 



(2) - u{k - u) 



[x-uixr 



(27) 



for S e Z>o. 

Proof: Markov inequality 

Pr[Cs{G) > 1] < E[G5(G)] (28) 

provides an lower bound on Pr[X{G) > S]: 

Pr[X{G)>6] = l-Pr[G5(G) > 1] 

> 1-E[G5(G)]. (29) 

■ 

IV. Numerical result 

In order to evaluate the tightness of the lower bound shown 
in Theorem |2] we made the following computer experiments. 
In an experiment, we generated 10*-instances of undirected 
graphs from the random graph ensemble defined in the SecjlL] 
|b] We assumed that q = 5 and = 0.1, /i(2) = 0.2, n{3) = 
0.4, /i(4) = 0.2, ^(5) 0.1; namely, 



(30) 



(25) 



The minimum cut capacity for each instance was computed 
by using the Ford-FuUcerson algorithm [15]. 

Figure [T| presents the accumulate distribution of minimum 
cut capacity when the number of vertices and edges are 
k = 100 and n = 400, respectively. The lower curve represents 
the lower bound presented in Theorem |2] and the upper curve 
is approximate values Pr[X{G) > S] obtained from 10"*- 
randomly generated instances. We can observe that two curves 
shows reasonable agreement in the range 1 < S < A. 

Figure |2]deals with a denser graph ensemble compared with 
that used in Fig. [T] In this case, two curves are very close in 
the range 1 < 5 < 15. Compared with Fig. [T] we can see that 
the proposed lower bound becomes tighter for a denser graph 
ensemble. 

From these these experimental results, it can be said that 
the proposed lower bound captures the accumulate distribution 
Pr[X{G) > 6] of the min-cut capacity of the random graph 
ensemble fairly well. 
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Fig. 1. Sparse case: accumulate distribution of tlie minimum cut capacity 
Pr[\{G) > S] (k = 100, n = 400): estimation by simulations and lower 
bounds 
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Fig. 2. Dense case: accumulate distribution of the minimum cut capacity 
Pr[X{G) > <5] (fc = 100, n = 700) estimation by simulations and lower 
bounds 

V. Conclusion 

In this paper, a lower bound on the accumulate distribution 
of the minimum cut capacity for a random graph ensemble 
is presented. From the compute experiments, it is observed 
that the lower bound reflects actual statistical behavior of the 
minimum cut capacity. The bound and the proof technique 
presented in the paper would deepen our understanding on 
typical behaviors of the minimum-cut capacity. 

The proof technique used here has close relationship to 
the analysis for the average weight distributions of LDGM 
codes with column weight 2. The one-to-one correspondence 
between a cut space and the row space of an incidence matrix 
implies that the minimum cut capacity is an analog of the 
minimum distance of the binary linear code defined by an 
incidence matrix. The analysis presented here has similarity 
to the typical minimum distance analysis of LDPC code 



ensembles [17]. 

An advantage of the proposed technique is its applicability 
for a graph ensemble with finite number of vertices and 
edges. Most related studies deal with asymptotic behaviors 
and cannot directly be applied to a finite size graph ensemble. 
Of course, it would be interesting to investigate the asymptotic 
behavior of the proposed lower bound when n and k approach 
infinity while maintaining the relationship n — f{k) (/ is a 
real-valued function, e.g., n = (3k^). 

The second advantage of the proposed technique is exten- 
sibility. In this paper, we discussed a simple graph ensemble, 
which is closely related to the Erdos-Renyi random graph 
ensemble [4]. The analysis for deriving the average cut-set 
weight distribution is approximately equivalent to the analysis 
of the average weight distribution of an LDGM code ensemble 
[6] or of the average coset weight distribution of an LDPC 
code ensemble [8]. Extension to other graph ensembles, such 
as regular or irregular bipartite graph ensembles, may be 
straightforward. 
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